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As an explanation for the acceleration of the universe, we propose dark 
matter with self-interactions characterized by a negative pressure; there is no 
vacuum energy whatsoever in this Cardassian model. These self-interactions 
may arise due to a long-range “fifth force” which grows with the distance 
between particles. We use the ordinary Friedmann equation, = SvrGp/S, 
and take the energy density to be the sum of two terms: p = pM + PK- 
Here pM = rnuM is the ordinary mass density and pK = Pk{pm) is a new 
interaction term which depends only on the matter density. For example, 
in the original version of the Cardassian model, px = bp'^ with n < 2/3; 
other examples are studied as well. We use the ordinary four-dimensional 
Einstein’s equations = SnGTf^iy and assume a perfect fluid form for the 
energy momentum tensor. Given this ansatz, we can compute the accompa¬ 
nying pressure, and use the conservation of the energy-momentum tensor to 
obtain modified forms of the Euler’s equation, Poisson equation, and continu¬ 
ity equation. With this fully relativistic description, one will then be able to 
compute growth of density perturbations, effects on the Cosmic Background 
Radiation, and other effects with an eye to observational tests of the model. 
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I. INTRODUCTION 


Observations of Type lA Supernovae [1,2], as well as concordance with other observations, 
including the microwave background [3] and galaxy power spectra [4], indicate that the 
universe is flat and accelerating. Many authors have explored possible explanations for the 
acceleration: a cosmological constant, a decaying vacuum energy [5,6], quintessence [7-9], 
and gravitational leakage into extra dimensions [10]. 

Freese and Lewis [11] proposed Cardassian expansion as an explanation for acceleration 
which invokes no vacuum energy whatsoever. In this model the universe is flat and accel¬ 
erating, and yet consists only of matter and radiation. They proposed additional terms in 
the Friedmann equation, which becomes = g{pM), where pm contains only matter and 
radiation (no vacuum), FT = d/a is the Hubble constant (as a function of time), G = 
is Newton’s universal gravitation constant, and a is the scale factor of the universe. The 
function q^Pm) returns to the usual SnpM /during the early history of the universe, 
but takes a different form that drives an accelerated expansion after a redshift z ^ 1. Such 
modifications to the Friedmann equation may arise, e.g., as a consequence of our observable 
universe living as a 3-dimensional brane in a higher dimensional universe [12]. 

In this paper we propose an alternative origin for the same behavior. This origin is a “hfth 
force” in ordinary four spacetime dimensions. We use the standard Friedmann equation 

= 87rGp/3, (1) 

but allow the dark matter to have self-interactions that contribute a negative pressure: these 
self-interactions drive accelerated expansion, and may arise due to a long-range “hfth force.” 
We speculate on a form of the force between particles that may be responsible for such an 
interaction: a conhning force. 

We take the energy density p on the right hand side of Eq.(l) to be the sum of two terms: 
the ordinary contributions from matter and radiation plus a new interaction term. During 
the matter dominated era, we take 
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P — Pm + PKi 


( 2 ) 


where pm = "nrnM (mass m times number density um of matter) is the ordinary mass density 
and pk is a new interaction term which is a function only of the mass density, px = Pk{pm)- 
We take the mass density to scale as usual with the redshift, 

Pm = Pm,o{o-/o-o) ^ (3) 


where subscript 0 refers to today. 

For example, consider the original Cardassian model proposed in [11]: 

Stt 


= 


3mli 


Pm + Bplj with n < 2/3, 


( 4 ) 


or, equivalently. 


= 


SttG 


Pm 


^ ( dCard I 

\ Pm ) 


1—n 


( 5 ) 


where pcard is the matter density at which the two terms are equal. The second term in 
square brackets is negligible initially (when pm 3> Pcard), and only comes to dominate at 
redshift 2 : ~ 1 (once pm ~ PCard)- Once it dominates, it causes the universe to accelerate: 
the scale factor grows as a ^2/3n ^ ^ 2/3 so that a > 0. 

There are two possible interpretations of Eq.(5). Previously, [11] interpreted it as a 
modihed Friedmann equation which arose from the physics of extra dimensions [12]. Here, 
and in the longer paper [13], we treat this equation as an ordinary Friedmann equation in 
which the second term describes self-interactions of the dark matter particles) with 


Pk — 


3mlj 

“sr®''"' 


( 6 ) 


In a ‘generalized Cardassian model’ other functions qIpm) of the matter density on 
the right hand side of the Friedmann equation can also drive an accelerated expansion in 


^This interpretation may have a four-dimensional origin or serve as an effective description of 
higher dimensional physics. 
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the recent past of the universe without affecting its early history [14,13]. Several of these 
alternative functions will be discussed below [see Eqs.(23) and (26)]. 

In the model we discuss here, we use the ordinary four-dimensional Einstein’s equations 

= SttGT^,. (7) 

We take the energy momentum tensor to be made only of matter and radiation, with the 
perfect fluid form 

^ ^ ^ p)u^^u’^, (8) 

where p is the total energy density given in Eq.(2), p = pm + Pk is the accompanying 
pressure, and is the fluid four-velocity. The Cardassian contribution to the pressure can 
be computed in our class of models and turns out to be a negative pressure, px < 0. This 
negative pressure is responsible for the universe’s acceleration. We can End models with any 
negative equation of state wk = Pk/Pk < 0, both constant and time-dependent^, including 
wk < —1; we also find models in which the pressure depends inversely on the energy density, 
p ~ -1/pE 

We use the conservation of the energy-momentum tensor to obtain modified forms of 
the Euler’s equation, Poisson equation, and continuity equation. With this fully relativis¬ 
tic description, one can then compute growth of density perturbations, effects on Cosmic 
Background Radiation and cluster abundance, and other effects with an eye to observational 
tests of the model. An initial study of density perturbations, performed in the Newtonian 
limit, can be found in our longer paper [13]. 

The Cardassian model also has the attractive feature that matter alone is sufficient to 
provide a flat geometry. The geometry of Eq.(l) is flat, as required by observations of the 

^We remark that in our class of models, wk < — 1 is not connected to a violation of the weak 
energy condition. Since px and px are not independent of pM, the weak energy condition is a 
condition on the total energy-momentum tensor: {pM + Px)/{PM + Pk) > ~1- 
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microwave background [3]. The numerical value of the critical mass density for which the 
universe is flat can be modihed; we take the value of the critical mass density to be 0.3 of the 
usual value. Hence the matter density can have exactly this new critical value and satisfy 
all the observational constraints, e.g., from the baryon cluster fraction and the galaxy power 
spectrnm. Throughont this paper, we will assnme that there is no vacuum energy at all. 
We do not solve the cosmological constant problem; we simply set it to zero. 

We speculate on the possible origin of an interaction energy with negative pressure in 
Sect. II, present a general fluid formulation in Sect. Ill, and then give specihc examples in 
Sect. IV. We conclude in Section V. 

II. ORIGIN OF INTERACTION ENERGY WITH NEGATIVE PRESSURE 

Here we speculate on a possible origin for an interaction energy with a negative pressure. 
Dark matter particles may be subject to a new interparticle force which is long-range and 
conhning. This force may be of gravitational origin or maybe a hfth force. It is the conhne- 
ment property that gives rise to a negative pressure (see [13]). That a conhning force can 
give rise to an effective negative pressnre is well-known in particle physics, where the MIT 
bag model is just such an effective description of quark conhnement. 

To be somewhat quantitative, let us consider a simple example of a power law interpar¬ 
ticle potential 


Uy = A,%, (9) 

where a > 0, is the distance between particles and A is a normalization constant. The 
corresponding force is a power law, F{r) oc The total new interaction energy of a 

system of N particles occupying a volnme of radius R will be Unew — AN'^R^, to within a 
nnmerical factor of order 1 dependent on the geometry. The total gravitational potential 
energy of the same system is, also within a factor of order nnity, Hgrav — where M 

is the total mass of the system. To play a cosmological role at the present time, the new 
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energy must be of the same order of the gravitational energy when R ~ Rh-, the current size 
of the horizon. Imposing that ~ f/grav at i? ~ gives us the normalization A = 
where m = M/N is the mass of a single particle. The magnitude of the new force on galactic 
scales is negligible compared to the gravitational force. This Newtonian formulation must 
be modihed at large distances because of the hnite speed of light and issues of causality. 
Notice that the Cardassian index n in the model is connected to the exponent a in the 
conhning force law through a = 3(1 — n). 


III. BASIC EQUATIONS 

A. Perfect Fluid 

As discussed in the Introduction, we use the ordinary four-dimensional Einstein’s equa¬ 
tions = STiGTniy where the energy-momentum tensor is made only of matter and radia¬ 
tion and has the perfect fluid form, + {p + p)u^u^. Here the total energy density 

p for matter includes not only the mass density pm (mass times number density) but also 
any interactions as in Eq.(2). 


B. Conservation laws and Evaluation of Pressnre 


The Bianchi identities guarantee the conservation of energy and momentum. 


= 0 . ( 10 ) 

In a comoving frame, energy-momentum conservation gives the (fully relativistic) energy 
conservation and Euler equations [15,16] 


P= (P + P), 


u 




p + p’ 


( 11 ) 

( 12 ) 
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where the dot denotes a derivative with respect to comoving time and the tensor = 
diMu — projects onto comoving hypersurfaces. We impose in addition that mass (or 
equivalently particle number) is conserved, 


= 0. (13) 

This will give us the usual dependence of the mass density on the scale factor of the universe, 
Eq.(3). 

From Eqs.(l) and (11), or equivalently from adiabaticity, we can hud the pressure due 
to the new interactions. 


Pk — Pm 



— Pk- 


(14) 


C. Newtonian limit 

Now we obtain the basic equations in the Newtonian limit. In Minkowski space, we 
write = 7 ( 1 ,n) and the metric 7 “^ = diag(— 1 , 1 , 1 , 1 ), where 7 = 1/yl — ^ and v is the 
3-dimensional fluid velocity. Then, = diag(—p, p,p, p). We consider a weak static held 

produced by a nonrelativistic mass density. Poisson’s equation becomes 

= 47rG(p-|-3p). (15) 

From = 0 we hud the continuity and Euler’s equations 

h^ + V.[(p + p)wl = 0, (16) 

^ + (17) 

at P + P 

Notice that we do not assume p p in the right hand side of Eqs. (15-17). 

With the additional constraint of particle number conservation, we have the ordinary 
continuity equation for matter, 

^ + V • (pmv) = 0. (18) 

Hence our basic nonrelativistic equations are Eqs. (15 - 18). 
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IV. THREE EXAMPLES OF CARDASSIAN MODELS 


A. Original Cardassian Model 


In the original Cardassian model of Ref. [11], 

( \ l-n 

PM ) 

with n < 2/3 as in Eqs.(4) and (5). The pressure associated with this model in the fluid 
approach follows from Eq. (14) as 


Pk = -(1 - n)pK- 


( 20 ) 


This model therefore has a constant negative wk = Pk/Pk = ~ Any observational 

test of Cardassian that depends only on a{t) thus has similar behavior to quintessence; 
however, tests that rely, e.g. on modified Poisson’s equations will have different observational 
consequences. 

One can now obtain the basic equations in the Newtonian limit for the p^ Cardassian 
model by substituting px and px of Eqs.(19) and (20) into Eqs.(15 - 18). In particular, the 
modihed Poisson’s equation becomes 


VV = 47rG 


Pm 


(2 - 3n) ( 


PCard 

V Pm , 


1—n 


Pm 


( 21 ) 


The modihed Euler’s equation of Eq.(17) gives rise to a new force 


dv 

dt 


Pm V Pm J Pm 


( 22 ) 


On galactic scales, this force destroys hat rotation curves (velocities tend to increase as one 
goes out to large radii in an unacceptable way). This problem for rotation curves persists in 
a relativistic generalization of this argument. The huid p” Cardassian case must therefore 
be thought of as an ehective model that applies only on cosmological scales. The examples 
in the following two subsections present alternatives to address the issues on galactic scales. 



B. Polytropic Cardassian 


Another class of models has 


P 


^internal “1“ PCard 


1 + 


PM y 

PCard / 


1 

9 


(23) 


with g 7 ^ 0. This model can be used on all scales (see below), but it does not quite £t 
the criteria of “Cardassian” since it does have a vacuum term. At late times in the future 
of the universe, when pM Pcard, this model becomes cosmological constant dominated 
with A = pcard- Eq.(23) is very similar to a model that was derived earlier [10] due to 
gravitational leakage into extra dimensions. 

Here, the pressure is 


P — Pm — PCard 


1 + 



(24) 


When the ordinary pressure pm can be neglected, this model obeys a polytropic equation of 
state 


P PCard 




(25) 


with negative pressure and negative polytropic index N = —1/q. For g > 1, the pressure 
varies inversely with the energy density (to the power g — 1). 

We must make sure that at the scales of galaxies and clusters the Cardassian pressure can 
be neglected compared to the ordinary pressure. At large matter densities, the Cardassian 
pressure is \pk\ — Pm^P ca^d- ^ galaxy or cluster with velocity dispersion a, the ordinary 
pressure is pm — Pm<^‘^- We want \pk\/pm — (Pcard/PM)^/^'^ 1- Taking a ~ 300 km/s 

and assuming px is unimportant out to ~ 100 kpc where pm ~ lO^Pcard, this condition 
amounts to g > 3. 

The models that match the rotation curves have pressure that scales inversely with the 
energy density. The condition g ^ 3 requires that p oc —^ where r = g — 1 ^ 2. This model 
is similar to the Chaplygin gas [18] which has r < 1; note that we differ from the Chaplygin 
gas in the value of the exponent that we find to be required. 
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Fits to supernova data [17] show that g ^ 2. Hence we hnd difficulties in making poly¬ 
tropic cardassian models, or generalized Chaplygin gas models, work at both cosmological 
and galactic scales. 


C. Modified polytropic Cardassian 


Let us return now to the original Cardassian proposal in which there is no vacuum energy 
whatsoever. A Cardassian model that can be used on all scales is 


P Pinternal A PM 


^ / PCard \ 

\ PM ) 


q{l-n) 


(26) 


For q = 1 this reduces to the original p” Cardassian model. The pressure follows as 


P = Pm - (1 - n)p 


'M 


1 + 


f PCard ^ 


/ PCard \ 
\ PM ) 


(27) 


\ Pm ) 

This model is interesting because the two parameters n and q are important on different 
scales. The parameter n sets the current value of tc ~ — (1 — n), and so can be chosen 
to fit the supernova data, while the parameter q governs the suppression of the Cardassian 
pressure at high densities, and can therefore be chosen not to interfere with galactic rotation 
curves and cluster dynamics. On galactic scales, the pressure again depends inversely on the 
energy density, but on large scales of the universe it depends linearly. Concrete comparisons 
with data will be presented in [17]. 


V. CONCLUSIONS 

An interpretation of Cardassian expansion as an interacting dark matter fluid with neg¬ 
ative pressure is developed. The Cardassian term on the right hand side of the Friedmann 
equation (and of Einstein’s equations) is interpreted as an interaction term. The total en¬ 
ergy density contains not only the matter density (mass times number density) but also 
interaction terms. These interaction terms give rise to an effective ngative pressure which 
drives cosmological acceleration. These interactions may be due to interacting dark matter. 
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e.g. with a long-range confining force or a fifth force between particles. Alternatively, such 
interactions may be an effective description of higher dimensional physics. 

A fully relativistic fluid model of Cardassian expansion has been developed, in which 
energy, momentum, and particle number are conserved, and the modified Poisson’s and 
Euler’s equations have been derived. In [13], a preliminary study of density fluctuations 
in the early universe has also been presented. There we developed a Newtonian theory of 
perturbations, but discovered curious gauge ambiguities in the relativistic theory that must 
be resolved in a future study. 

One of our goals is to allow predictions of various observables that will serve as tests 
of the model. The Cardassian model will have unique predictions, particularly due to the 
modified Poisson’s equation. For example, one can now calculate the effect on the Integrated 
Sachs Wolfe component in the Cosmic Microwave Background, as well as the effect on 
cluster abundances at different redshifts. These predictions can then be tested against 
measurements of these quantities. Comparison with existing and upcoming supernova data 
is being studied in another paper [17]. 

ACKNOWLEDGMENTS 

We thank T. Baltz, J. Frieman, and Y. Wang for useful conversations. K.F. acknowl¬ 
edges support from the DOE via the University of Michigan, and thanks the Aspen Center 
for Physics for hospitality. P.G. thanks the Michigan Center for Theoretical Physics for 
hospitality. This research was supported in part by the National Science Foundation under 
Grant No. PHY99-07949 at the Kavli Institute for Theoretical Physics, Santa Barbara, CA. 


[1] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 
(1999) [astro-ph/9812133]. 


11 




[2] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) [astro- 
ph/9805201]. 

[3] C.B. Netterfield et al, astro-ph/0104460; R. Stompor et al, astro-ph/1015062; N.W. Halverson 
et al, astro-ph/0104489; C. Pryke et a/Astrophys. J. 568, 46 (2002) 

[4] L. Verde et al. [2dF Survey] MNRAS 335, 432 (2002) 

[5] K. Freese, F.C. Adams, J.A. Frieman, and E. Mottola, Nucl. Phys. B287, 797 (1987). 

[6] J. Frieman, C. Hill, A. Stebbins, and I. Waga, Phys. Rev. Lett. 75, 2077 (1995). 

[7] L. Wang and P. Steinhardt, Astrophys. J. 508, 483 (1998). 

[8] R. Caldwell, R. Dave, and P. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998) 

[9] G. Huey, L. Wang, R. Dave, R. Caldwell, and P. Steinhardt, Phys. Rev. D59, 063005 (1999) 

[10] C. Deffayet, G. Dvali, and G. Gabadadze, Phys. Rev. D65, 044023 (2002). 

[11] K. Freese and M. Lewis, Phys. Lett. B540, 1 (2002) [astro-ph/0201229]. 

[12] D. Chung and K. Freese, Phys. Rev. D 62, 063513 (2000). 

[13] P. Gondolo and K. Freese, hep-ph/0209322 

[14] K. Freese, hep-ph/0208264. 

[15] S. Hawking, Ap. J. 145, 544 (1966). 

[16] D. Lyth and E. Stewart, Ap. J. 361, 343 (1990). 

[17] Y. Wang, K. Ereese, J. Erieman, P. Gondolo, and M. Lewis, in preparation. 

[18] A. Kamenshchik, U. Moschella, V. Pasquier, Phys. Lett. B551, 265 (2001); M.C. Bento, O. 
Bertolami, and A.A. Sen, astro-ph/0210468; N. Bilie, G.B. Tupper, and R.D. Viollier, astro- 
ph/0207423; V. Gorini, A. Kamenshchik, and U. Moschella, astro-ph/0209395. 


12 



